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SUMMARY 


The aim of this research is to recast the indicial response approach in a form appropriate to the study 
of vortex-induced oscillations phenomena. We demonstrate that an appropriate form for the indicial response 
of the velocity field may be derived directly from the Navier-Stokes equations. For convenience, the study is 
divided into three parts. In Part I, we demonstrate, on the basis of the Navier-Stokes equations, how a form of 
the velocity response to an arbitrary motion may be determined. To establish its connection with our previous 
work, the new approach is applied first to the simple situation wherein the indicial response has a time-invariant 
equilibrium state. Results for the aerodynamic response to an arbitrary motion are shown to confirm the form 
that we have obtained previously. 


INTRODUCTION 


The motivation behind this research originates with a recent review article by Tobak et al., which is 
an exposition of their approach to the modeling of the aerodynamic contribution to the inertial equations of 
motion of a maneuvering aircraft (ref. 1). The approach features the use of nonlinear indicial responses and 
generalized superposition integrals. Recognizing its natural connection with ideas from bifurcation theory, the 
authors of reference 1 extended their modeling approach to accommodate the potential occurrence of bifur- 
cations involving time-dependent equilibrium states. Of particular interest is the first of these, called a “Hopf 
bifurcation,” wherein a formerly stable time-invariant equilibrium state is replaced by a time-varying periodic 
equilibrium state. In the implementation, however, Tobak et al. discovered that an important additional issue 
arises along with a Hopf bifurcation: the necessity of specifying phase as well as amplitude and frequency 
to completely determine the periodic equilibrium state. Although the amplitude and frequency of the peri- 
odic state are determined entirely by flow conditions that are independent of an origin in time, specifying the 
phase requires additional information involving an initial condition. Through specification of the phase, the 
equilibrium state is required to acknowledge an origin in time. 

In aerodynamic applications, the physical origin of a large-scale periodic equilibrium state often is 
the onset of vortex-shedding. Of the many examples, Tobak et al. (ref. 1) cited the occurrence of stall on an 
airfoil when the angle of attack exceeds a critical value, and the wake of the flow past a cylinder when the 
Reynolds number exceeds 50. The latter example, vortex-shedding from a cylinder or, more generally, a bluff 
body, is itself a subject that has attracted the attention of researchers for many decades. When the body is 
able to interact with the vortex-shedding through, for example, elastic mountings, the consequences are varied 
and potentially dramatic. This subject, which has come to be known as ‘‘vortex-induced oscillations,” has 
also inspired a large literature, reflecting the wide range of practical circumstances in which structures may be 
endangered by undergoing such oscillations. Excellent surveys of the field have been published by Sarpkaya 
(ref. 2) and more recently by Bearman (ref. 3). Both authors noted that in spite of the long and intense interest, 
an adequate rationally based mathematical model of the phenomenon still does not exist. 

In a first attempt to fulfill the need for a suitable mathematical model of the aerodynamic contribution 
to the equation governing the motion of an elastically mounted cylinder immersed in a uniform oncoming 
stream, the authors of reference 1 proposed a simple model based on one of the alternative approaches studied 
therein. However, the model does not incorporate the memory effects that the specification of phase would 
have required. A careful numerical study of the resulting equation of motion, carried out in reference 4 and by 
us, suggests that the model proposed is unable to predict the occurrence of the self-sustained oscillations that 
experiments have revealed. 
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To that end, this study has two purposes: (1) to show how the deficiency of the model proposed in 
reference 1 may be remedied to include the effects of memory; and (2) to apply the results of the analysis 
to the particular case of flow past a cylinder that is in periodically forced transverse motion. As it turns out, 
revising the reference 1 model to include memory effects must be done at the level of the velocity field itself. 
We demonstrate that a form for the indicial response of the velocity field may be derived directly from the 
equations governing the fluid motion. These are taken to be the Navier-Stokes equations for time-dependent 
incompressible flow. 

For convenience, the study will be divided into three parts. In Part I, we demonstrate on the basis 
of the Navier-Stokes equations how an indicial response for the velocity field may be formulated, and then 
used within a generalized superposition integral to determine the form of the velocity response to an arbitrary 
motion. To establish its connection with our previous work, the new approach is applied first to the simple 
situation wherein the indicial response has a time-invariant equilibrium state. Results for the aerodynamic 
response to an arbitrary motion are shown to confirm the form that we have obtained previously by a variety 
of approaches (ref. 1). In Part II, the analysis is then directed to the new situation where the equilibrium state 
of the indicial response is periodic in time. A method is derived for specifying its phase, and, by allowing a 
correct incorporation of memory effects, this proves sufficient to remedy the deficiency of our previous model. 

Part III is the object of the second purpose, namely, to apply the results of the analysis to the particular 
case of flow past a cylinder that is in periodically forced transverse motion. We show that our approach 
captures the distinctive features of vortex-induced oscillations that have been revealed by the results of careful 
experiments on this motion (reported in refs. 2, 3). Finally, we show how the indicial response approach may 
be reconciled with a currently popular approach based on the use of amplitude equations. 

GENERALIZED SUPERPOSITION INTEGRAL BASED ON 
INDICIAL RESPONSE APPROACH 


We derive the indicial response approach with emphasis on the physical postulates involved. 


For generality, consider an aircraft that has started from rest in the distant past with fixed axial velocity 
Uo and zero vertical velocity (cf. fig. 1). Its motion is referred to an X, Y coordinate system that is fixed in 
space. It passes through the origin at the arbitrarily chosen initial instant £ = 0 , maintaining the constant axial 
velocity Uo and simultaneously translating vertically, with the vertical velocity v c at the center of gravity being 
an arbitrary function of time £. The angle of attack cv is defined as the angle between the resultant velocity 
vector and the aircraft’s longitudinal axis: 


a = tan 


-\ 



( 1 ) 


Let us note that we specify a constant axial velocity Uo to be in accord with normal operating conditions in 
wind- or water-tunnel experiments, wherein the uniform oncoming flow, normally held at constant velocity, 
would supply the corresponding value of Uo « To form the indicial response, we need to consider two motions. 
In the first one, the aircraft undergoes a variation of angle of attack c*(£) from time zero to a time £ = r (cf. 
fig. 2). Subsequent to time r, the angle of attack is held constant at a(r) . The first motion history therefore is 
designated a r (£): 


c* r (£) = 


<*(£) 

Qf( T) 


if 0 < £ < r 
if i > r 


( 2 ) 
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Figure 1. Maneuver referred to space-fixed ( X , Y) and moving ( x , y ) coordinates, the latter attached to the 
airplane. 


«t(£) = «(€); 
t\ =o(x);6>x 


«i(9 = «(*); 

/\ = a(t)+Aa; ^ S=t 





The angle-of-attack history a T (0 is said to belong to a family of motions, all of which are generated 
by the “natural” motion history or(0, and all of which have the characteristic of remaining constant beyond 
the instant designated by the subscript, here r. Accordingly, the angle-of-attack history ov+atCO belongs 
to the same family, being the extension of the history ar(0 to the instant r + A r, beyond which a remains 
constant at or( r + A r) . The second motion needed to form the indicial response, however, does not belong to 
this family. In the second motion, the aircraft undergoes the same angle-of-attack history ar(0 up to time r. 
Subsequent to r, the angle of attack again is held constant, but is given an incremental step change A a over 
its previous value of a( r) . The second motion, designated ar*(0 , is represented as 

o£(O = ar(O + Aa..ff T (0 (3) 

where H T ( 0 is the Heaviside function: 


H t (0 = 


0 

1 


if 0 < £ < r 

if i>T 


(4) 


Let us define by u(x, t, t)(x : spatial coordinates, t: time of observation subsequent to r) the velocity field in 
the neighborhood of the aircraft. 


The indicial response approach is established by making use of the following postulates. 

POSTULATE 1: Corresponding to a motion history belonging to the family of motion histories a r (£) 
with r 6 [0 , £]> there exists a velocity field u[af, a r (£); t, r] well defined and supposed known by some 
means. 


POSTULATE 2: To each member of the family of motion histories a£(0 with r 6 [ 0 , t] , there exists 
a velocity response u[ x, oi*(Q ; t, r] . For any observation time t > r + A r, this velocity response reproduces 
the velocity response corresponding to the motion history of r+ Ar(0 within a negligible error of i?( A a ) 2 : 

u[x,a*(0] =u[£,a T (0 + Aa.tf T (0] = u[f,cw T (0] + tf(Aa) 2 (5) 

The increments A a and A r are connected by the rate-of-change of the angle of attack a( t) = ^ , which is 
presumed to be defined for each value of r E [ 0 , t ] . 

POSTULATE 3: For every value of r 6 [ 0 , t] , there exists an indicial response u a [ x, a (() ; t, r] 
defined as 

Uo,[i,a( 0 ; t,r] = lim Aa -> 07 — {u[?,ar (0 + Aa.H T (0\ t,r] - u[ x, o: r (0 ; t, r] } (6) 

Aar 


Given a motion history at(£ : £ E [ — 00 , i]), it is possible to choose an initial time tq and partition 
the time interval [ tq , t ] into [ to , ti , . . . , r n ] such that 


T n = t , tiAt = t , Ti - 7?_1 = At (i= 


One has the following relations: 

ui[ar 7 i ,(0 + Aar • ff TO (0] - tfotor^CO] 

«2[cv7i(0 + Aar • H n (0] - t?t [ afr, (01 


u«[ a T ^ (0 ] + A or • H r ^ (() ] - u^i [ (0 1 


tfa[Qr 7O (0]Aar+ i?(Aar) 2 
Uatar^ (0]Aar + tf(Aar) 2 

41^1 (0]Aar+ i5(Aar) 2 J 


(7) 


( 8 ) 
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Note in (8) that the additional dependencies of u, on x, t , t,( » = 1 n) have been omitted for brevity. 

By summing up the relations (8) and using Postulate 2, one gets 

n-1 

u n [a TB (O]-t? 0 [a TO (O] = X)^^(O]Aa+t?(Acv) 2 (9) 

i=0 

or, in the limit as A r — > 0 , 

u[x,a ( (0'.M] = tfl i, 0^(0^170] + J dr u a [i, a r (0; t, r]^ + i?(,Aa) 2 (10) 

The relation (10) constitutes our main result that will be used in subsequent sections. Notice that it relies on 
three physical postulates which have to be satisfied. 

CONSTRUCTION OF VELOCITY RESPONSE FOR 
MOTION HISTORY: a T (0 

In order to derive the aerodynamic response on the basis of relation (10), we need to establish forms 
for the velocity responses to the motions a T (() and a*(Q over the time interval t — r > 0 . In this and the 
following section, we shall demonstrate that the forms can be derived analytically from direct consideration 
of the Navier-Stokes equations. 

For simplicity, we neglect compressibility effects and assume that the fluid motion is governed by the 
Navier-Stokes equations for an incompressible fluid. In a coordinate system attached to the body, the Navier- 
Stokes equations have an additional term v c (where v c is the vertical velocity at the mass center and the dot 
denotes a derivative with respect to time) to account for the acceleration of the coordinate system relative to 
inertial space: 

c\ -* 

^ + (u.V)u +Vp-vV 2 u = -i? c (11) 

at 

Here, u is the velocity field, p is the pressure (normalized with respect to constant density p), and v is the 
kinematic viscosity. We wish to derive a form for the velocity response to the motion a T (0 over the time 
interval t — t > 0. Since a T (0 = a(r) = constant for t — t > 0, boundary conditions determining the 
velocity response over the interval t - t > 0 are perfectly steady. Consequently, we expect the velocity 
response to attain an equilibrium state as the interval becomes large, that is, as t — r — ► oo. 

The response of the velocity field can be decomposed into an equilibrium state u equ ii and a transient 
component utrans that decays as time increases. The principal condition that we imp.se in Part I is that the 
equilibrium state u eqil u be time-invariant. Thus, as time is referred to r (the instant after which the angle-of- 
attack is to be kept constant), we have 


u( X , t — T > 0) 

— Uegut’fC?) + Utrans T > 0) 1 

(12) 

/-—S 

C“+* 

1 

H 

V 

O 

“ Pe^tiiK + Ptrcms T>0) J 

To simplify the notation, we define 


t+ = t — T > 0 

(13) 


Let us assume that the equilibrium state (u e?ul /( i) , p eq u<i( x) ) is available, for example, from a solution of the 
time-independent ( = v c = 0) version of equation (11). 
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Substituting ( 1 2) into (11) and noting that tf c = 0 at t + = 0( +) , we get the following for the equations 
governing u, rtms : 


Utroru + ( Uejuti-^) U( ronj + ( Ut rans . V) U e ^ u ,J + Vp^rotis U trails + ( Uj rans .V) Ui rans — 


V.tftrons - 



with u trana = 0 on appropriate boundaries d£2 . 


To find a form for ut ron s . we shall first consider a small disturbance u which satisfies the linearized 
version of equations (14) 

+(u effu ,i.V)u +(u.V)u e , ulJ + Vp-i/V 2 u = 0 1 

V.u = 0 j ^ V 

with u = 0 on dCl . Equations (15) constitute a linear eigenvalue problem. Since u equ a is independent of 
time, the eigensolutions 7 „ are also independent of time. It is possible to choose them such that (from ref. 5) 

u = u n = e x " u 7 „(i) , V.'fn = 0 , 7 n| da = 0 (16) 

From equations (15), the corresponding expression for p is 

p=e x ' u p n (x) (17) 


The eigensolutions 7 „( x) ( n = 1 , . . .) are associated with the eigenvalues \ n , the latter all having negative 
real parts in the physical situation under consideration, wherein the equilibrium state u e?u »i( x) is supposed to 
be stable to small perturbations. There is an equation adjoint to (15) having a set of eigensolutions such 
that (from ref. 5) 

{l[, li) = f v d,Vij(x).7i(x) = 8ij Vi,; 1 ng . 

(T^,7,) = f v dV 7j(x).ii(x) = 0 Vi,; J 

where the brackets ( ) denote the scalar product over space. The eigensolutions 7 „( x) span a complete 

functional space, and one can use this fact to construct a suitable solution form for ut rans (x,t+) . Returning 
now to the full nonlinear equations (14) governing u trons , let us assign u t Tana the form resulting from its 
projection onto the functional space of the 7 n : 

Utrons(f,t+) = E ( rf »( t +)'7n(l) + <*».(*+) 7, „(x)) (19) 

n 

where the barred terms denote complex conjugates (abbreviated c.c.). It will be convenient to break the pres- 
sure term pt ron * in (12) into two parts: one following the form of equation (19) and correponding to the linear 
part of the disturbance » ,uations, and the other (denoted p) corresponding to the nonlinear part of the distur- 
bance equations: 

PtTons(x,U ) = 52 (dn(U)p n (x) + d„(U)f n (x )) + p(l,t+) (20) 


Alter inserting equations (19) and (20) into (14) and eliminating the set of terms that identically satis- 
fies the linear form (15), we find that the coefficients dn(t+) satisfy the following equation: 

5 . Un y ^ ( dn ^n^n)7n ^P = ^ 1 ^ { ^n < ^m(7n-^)7»n+ dm ( Un(. r 1n-'^)lrn 

n n n ,m 

+^^(7n-V)7m + 4^(7„V)7 m ) (21) 
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Multiplying equation (21) by the vector 7 ^( 1 ), integrating over space, and using the properties of adjoint 
vectors defined by (18), we obtain 

dj(U) - Xjdj = -52 (7n-V)7m> + (7« V)f ro > 

n,m 

+dnd m (Tj,0 n .V)i m )+ ; = ',■■■ (22) 

Notice that the pressure gradient term Vp in equation (21) makes no appearance in equation (22). The prop- 
erties of 7 „ and 7 ^ (solenoidal and vanishing on < 9 £i ) ensure that the scalar products they form with gradient 
terms (such as Vp) will be identically zero. By defining appropriate coefficients A }nm , B }nm , C ;Tm » and D jnm , 
we rewrite ( 22 ) as: 

dj(U) - \j(a(r))dj(U) = 52 (Ajnm(o‘('r))d n (t + )d m (U) + B; wn (a(r))ci n (f + )3 Tn ( f+) 

n,m 

+ Cjrvm ( &C 1"))d n (t+) d m ( t+ ) 

+ D jwn {a(T))d n (U)d n (U)) j = 1,... (23) 

Equations (23) are transformed into integral equations: 
dj(U) = eM«(r))U dy(t+ = o) + 52 A ;nro (a(r)) f* ds e x ^ T)Uu ~ s) d n (s)d m (s) 

n,m ® 

+ 52 %»»(“( r » r dS e X,iaiTmU ~ S) dn(s)d m (s) 

n,m 

+ 52<W°< r » r ds ^ (aiT)Hu -* ) d n (s)d m (s) 

+ 52 £>;«*»(<*( r)) f* ds e^ ( “ (T))(t+ - s) 3 n ( S )d rn (s) j = 1 , . . . (24) 

n,m 


Solutions of the system of equation (24) can be obtained by using Picard’s method of successive 
approximations: 


(* + ) = 
^(M = 


e Mo <r))U djiu = Q) 

4 0) (i + ) + 52>W“(r)) f U dse^ t *~ a '>d^(s)d^Hs) 

n,m 


= 4 n) (^) + E^(^ T ))/ t+dseM ‘ + " a) 4 n) («)4 n) (5) 


where for brevity we have omitted writing the nonlinear terms in dndm , dnd m , d n dm- 


( 25 ) 
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Solutions dj(t + ) of equations (25) have the form of functional expansion series: 


dj(U) = d ( - 0 \u) + 'E,Aj rm (c l (T)) f' f d3oe^ u - 3o) S°Hso)d^Hso) 

+ y ]A im n(a(T)) f U dso H ds i e x ' (u - S0) \S n 0) (so).e Xm(a °- st) £ A mpq (a(T)) 
n,m J ° J ° l M 

x4 0) (3 1 )4 0) (5i) + e x -* (ao - i, > ) X;A w (a(T))4 0) (5,)d; 0) (3i).cC ) (so) 

P)9 

rt+ rso rs\ 

+ ...+ / dso / dsi / ds2 • • • (26) 

Jo Jo Jo 

The convergence of the series is insured by imposing, for instance, a Lipschitz condition on the integrands 
of the four integrals in the right-hand side member of equation (24). Calculation of the series lends itself to 
symbolic computation, as shown in reference 6. 

To complete the solution for the dj(t+)(j = 1 , ) at any level of approximation, we are required 

to specify their initial values at t+ = 0. The values of d ; (t+ = 0) must be determined from a match of the 
velocity fields on either side of £ = r: 


u(f,£ = r_) = u(£,£ = r+) (27) 

On the negative side of £ = r, in general the velocity field is the resultant of the entire history of the motion 
c*(0 up to the “present” time £ = r. Through specification of the dj at t+ = 0 in terms of this velocity 
field, the velocity response for t+ > 0 acknowledges its dependence on the past motion. The representation 
of this past motion in a suitable way is all that remains to be done to realize a complete characterization of 
the velocity response for > 0 . We shall defer a discussion of this important step until we have derived the 
velocity response to the motion cv*(£) and formed the indicial response. For the present, we simply designate 
the dependence of d/(t+) on the past motion as a functional 

dj(U) = dj[otU:te[0,T)) ,U] • (28) 

and write the velocity response to the motion a r (0 for t+ > 0 in the form 

uti.arCOU.r] = Uequii(S,a(T)) + (£„[<*(£ € [O.TDjt+l^nCx.aCr)) + c.c. (29) 


CONSTRUCTION OF VELOCITY RESPONSE TO 
MOTION HISTORY: a*(0 

In addition to constructing the velocity response to the motion a*(Q containing a step change in a 
at £ = t, we need to be assured that the derived response will satisfy Postulate 2. To that end, consider the 
interval t <(< t+ Arm which the velocity response changes continuously from u[ a T ( Q ] to til «r+Ar(0] 
as a result of the continuous change in a over the interval. We intend to show that if the continuous change 
in a is replaced by a step change, applied anywhere between r and r + A t, the resulting velocity response at 
any observation time t > r+ At will differ from the response to the continuous change in a, u[a T +Ar(0]. 
by terms of t?( A a) 2 . 
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In the equation of motion (11), for brevity let 

£[u] = (u.V)u + Vp-i/V 2 u (30) 


and let uj be the velocity response to the continuous change in cv for r < £ < r + A r. With u = uj and with 
the time variable t+ = t — r, integrate equation (1 1) from t+ = 0 to t+ = A r 


which gives 


Mr 

Jo dt+ 


tiidt+ 


uj(Ar) — ui(0) + 


r 

Jo 




v c dt+ 


C[u\]dt+ = — (v c (A r) — v c (0)) = -Au c 


(31) 

(32) 


Now let the continuous change in a over the interval r<^<r+Arbe replaced by a step change of the same 
net amplitude A a = o(r+Ar) — a( r) , applied at an arbitrary time r* anywhere between r and r + A r. Let 
U 2 be the corresponding velocity response over the interval 0 < t+ < At, and integrate equation (11) again 
(note that the initial state at t+ = 0 remains unaltered so that U 2 (0) = u\ (0)): 


[ ^7-U2 dt+ + f C[u 2 ]dt+ = -Av c [ 8(t+ -(T* — T))dt+ 

Jo at+ Jo Jo 


(33) 


which gives 

U2(At) — U2(0) + [ C[u2]dt+ = —Av c 

Jo 

Taking the difference between equations (32) and (34) yields 


U 2 (At) — ui(Ar) + f (£[u 2 ] -£[ui]) dJt + = 0 

Jo 


But 


Hence, in equation (35) 


ui = u[o(t(0] + terms of t?(Aa) 
U2 = u[<2t(0] + terms of i)(Aa) 


[ (£[u 2 ] - £[ui]) dt+ = ^ (£[u 2 ] - £[ tti ]) = tf(Aa) 

Jo a 


(34) 

(35) 

(36) 


(37) 


where the underline denotes an average over time. 


Therefore, U2 and ui differ by terms of i?(Aa) 2 at t+ = At. For t+ > At, note that ui = 
ui a r + A t ( 0 ] . whereas u 2 = u[a^*(0]. Since for all t+ > At, boundary conditions and equations of 
motion are identical for the two velocity responses, their solutions will differ only by terms of the order of 
the difference in initial conditions, which, as we have just seen, is of il(Aa) 2 . Thus, we are assured that the 
velocity response u[a£>(0] will satisfy Postulate 2. For simplicity, we put, henceforth, r* = t. It remains to 
develop a form for u[a*(0] analogous to that foru[o! r (0]. 


The derivation for tf[ a*( 0 ] is required to account for the effect of a step change in a at £ = t. This is 
reflected in the equation of motion by the presence of an impulsive forcing term arising from the time-derivative 
of the step change in v c . 


du - — i 

~r— + (u.V)u + Vp - uV u = -v c = -Av c 8(U) 

(7t+ 


(38) 
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where, to first order in A a, the magnitude of the step change in v c , denoted A v c , is 

A v c = Uosec 2 (a(T))Aa (39) 

The derivation parallels that for u[ a r (£)]. Again, the velocity response is decomposed into a time-invariant 
equilibrium state and a transient component, 

u( X, ) = U egui i( x, a(r) + Aa) + u trons ( x , U , a( r) + A a) ( 40) 

where now both components are to be evaluated at the new level of angle of attack, a( r) + A a. The equation 
for tftrons takes the same form as that of equation (14) with the exception that the right-hand side now contains 
the impulsive forcing term — A v c 8( t+ ) . 

As before, we try to form a solution by a suitable projection onto the functional space of the set of 
eigensolutions of the linearized homogenous equation 

Wtrons(l,U) = £n f ^n( U) 7»( l) + Mt+Hn(z) 

Ptrans ( I, t+) = £n (M *+)?'„(*) + h n (t + )p n (x) 

where the eigensolutions (7 n ,p n ) 316 now those appropriate to the new level of a, a(r) + A a. Multiplying 
the resulting equation (the counterpart of eqs. (21)) by the adjoint vector -f*( x, a(r) + A a), integrating over 
space and simplifying, we get the counterpart of equations (23): 

hj(.t+) — ) 1 (Ajnmh n h m + . . . + Dj rm fi n h m S J — —A a cj6(t+) J = 1 , . . . (42) 

n ,m 

where 

Aa c ; (ck(t) + Aa) = J ^ *(x,a(T ) + Aai).Av c dV (43) 

and the eigenvalues as well as the coefficients A ;nm D ;nm have the same definitions as before, albeit, 

now referred to a( r) + A a rather than to a( r) . 

Writing equation (42) as integral equations yields 
hj{U) = C } e x,t * — Aa cjH(U)e Xft * + J2 A in* ds e x ’ <u - s) h n (s)h m (s) + ... 

+ E Djm* r 4* e x ^~ a) h n ( s)h m (s) j = 1 , . . . (44) 

^ Jo 

which again can be solved by successive approximations, once the open constants Cj have been determined. 
To determine them, a condition on the velocity fields on either side of £ = r must be imposed. Here, it is 
no longer appropriate to require the continuity of the velocity fields across £ = r. Rather, we must admit the 
existence of a discontinuity in the velocity fields across £ = r, reflecting the step change in a that occurs there. 
Accordingly, we shall require that the velocity fields on either side of £ = r differ by an amount proportional 
to the step change in a at £ = r. The condition is 

trK(0.£ = T + J -u[c* t (£U= T_] = AaH(0)6(x,T) (45) 


+ p(f,t+) 


(41) 


10 



Using equations (12) and (19) for u[a T (Q,£ = r_], and equations (40) and (41) for u[ «*(£),£ = t+], we 
have in equation (45), 

(ue,ut l(x,a(T) + Ack) - Ueqml(x,a(T))) + E (^n(f+ = 0)%(x,a(T) + A a) + linin') 

n 

-E (<*»(*+ = 0 Hn(SMr)) + dnin) = Ac* H{U = 0) O(x.r) (46) 

n 

Multiplying equation (46) by the adjoint vector 7 ; ( x, a( r) + A a) and integrating over space, we get 

hj(t+ = 0) = (^(a + Aa),-(tt e ,„ii(o(+ Ac*) - u equ ii(a)) 

+ E (<*»(*+ = 0 ) 7 n(a) + + Ac* /f(f+ = 0 ) <h(x, r) } ( 47 ) 

n 

Comparing equation (47) with (44) evaluated at i+ = 0 , we have 

Cj = ( ij(a+ Aa),- (ueguiK ot+ A a) - u eju ,i(a)) + E (<*»(*+ = °)7n(c*) + d„7„) ) (48) 

n 

and 

( ij(x,a(T) + Aa),0(x,r) ) = -c ; (49) 

Equation (48) determines the form of Cj, and a comparison of the form of cj (cf. eq. (43)) with equation (49) 
reveals that A a O is simply A v c . The latter result implies that the velocity step imposed to the body, through 
the step change in angle of attack, is transferred instantaneously to the velocity field. The form of Cj can be 
simplified by retaining only terms to the first order in A a, since that is all we shall need to form the indicial 
response. After making extensive use of the properties of scalar products of adjoint vectors, we get 

C; = d ; (t + = 0) + Aa^d ; (t + = 0)+P^ + i?(Aa) 2 (50) 

with 

Vj = -{ij(x,a(T)), J^u[e* r (0,£ = r]) (51) 

Returning to the integral equation (44) for the hj(t+) with the values of Cj from equation (48) inserted, 
we get the following equation where only terms to the first order in A a are retained: 

hj(t + ) = e x ' ( “ (T))t+ dj{U = 0) + Ac*. {e x * a{r))u dj(U = 0)) + e x ' (a(r))i+ (£>, - H(U)cj) j 

+ E P ds{Aj nm e x ' (u - a) h n (s)h n (s) + ...+ Dj nm e x * u - S) h n (s)h m (s)} + Aa.E T ds 

x (A jnm e^ u - a) ) h n (s)h m (s) + ...+ ^ (. Dj nm e x * (u ~ 3) ) h n (s)hM j 

+ i?(Ac*) 2 (52) 

Anticipating the form of the result, we let 


hj(t+) = d ; (t+) + Ac*. ( — d ; (t+) + fj(U) ) + i?(Aa)' 


(53) 
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Inserting equation (53) in (52), we get linear integral equations for the /;(t+)0 = 1 , . . .) of the form 
fjiU) = (pj - CjH{U )) e X/t * + f e x * u -> (d„(s)f m (s) + d m (s)f„(s)) ds 

n,m 

+ /‘V' ( ‘*- a) (i(«)7*(«) + 3n.(*)7.(*))<fa (54) 

n,m 

where the eigenvalues A ; and the coefficients A ;nm D ;im are to be evaluated at a( r) . Equation (54) can 

be solved by successive approximations. Given the form (53) and assuming the f } determinable from (54), 
we are now able to write the form of tf[ a *( £) ] (cf. eq. (40)) as 

uTa*(0] = u egu n(x,a(r) + A a) + (dn(t+) + A aj^dn(U) + Aa/„(( t )j 7 „(f ,a(r) + A or) + c.c. 

(55) 

Again systematically retaining only terms to the first order in A a in (55), we have finally, 

*[<*<©] + E (/n( U ) 7 n( X, a( r) ) + / n 7 „) j ( 56 ) 

INDICIAL AND TOTAL RESPONSES OF THE VELOCITY FIELD 


«[«?(©) = t*IOf T (©]+ A 


- i 


Hie forms we have derived for u[ a T (£) ] in equation (29) and u[ a*(£) } in equation (55) satisfy 
Postulates 1 and 2. Postulate 3 presumes the existence of a limit as A a — ► 0 such that the indicial response 
can be defined as 


u Q [ i, a (0 ; L t] 


lim 

A a— 4 ) 


, Aa ) 


(57) 


We see from the way we have derived equation (56) that the limit (defined by the right-hand side of eq. (57)) 
will in fact exist, so that Postulate 3 can be satisfied. From equation (56), we have for the indicial response of 
the velocity field the following: 

u a = —u(x } ot(T);t,T) + 'S2f n (U)in(x,at(T)) + c.c. (58) 

da V 


The total response of the velocity field is obtained by summing the indicial response along the motion history 
at(0 according to equat : on (10) 


u( X, £) — u(x, 7{) ) + ^ d,T | dot ^ irons ) "t" ^ ^ /»( U ) 'Yn( x, a(_ t ) ) 


+ c.c. 


(59) 


that is, 


u( x,t) = u(x,Tt)) + u eguit ( x, a(t)) — U egu ,i( x,a(To)) 


/ \lLi ^52 dn ( i+ )'7n( ;? . Qf ( T )^ + ^,fn(.t + )i n (x,a(T)) + c.c.j (60) 


12 



In experiments, one usually maintains the angle of attack of the system constant for some sufficiently 
long time before letting it oscillate. Under these circumstances, the velocity field at £ = 7q is practically equal 
to its equilibrium value. Therefore, equation (60) is reduced to 


u(x,t) = u equil (x,a(t )) + J a(r) j + ^/ n (t+)7„(x,a(T)) + c.c.j 

(61) 

Equation (61) has the following form: 


r* da - r 

u(x,t) = u e quil(x,a(t)) + J dr— T [d, [a T (£) ,U = 0] ;t+] (62) 

with the functional dependence of T determined by the initial condition d } (t+ = 0)(; = 1 , . . .) , according to 
equations (25), (26), and (54). 


DETERMINATION OF THE INITIAL CONDITIONS dj(t + = 0)(; = 1 , . . .) 


According to equations (25), (26), (54), and (61), the indicial response of the velocity field is com- 
pletely known if the initial values d,( f+ = 0) (; = 1 , . . .) are determined. It is possible to derive an equation 
governing the behavior of d ; (f+ = 0) . Indeed, as stated previously, the value of the velocity field at £ = r is 
the resultant of the entire motion history up to the present time ( = t, a T (( : ( e] — oo, r[). According to 
relation (10), it is equal to 


u(x,£ = r) = u(x, ( = — oo) + f ds a(s) u a 

J — OO 

By using equations (29) and (61), equation (63) can be transformed as 
Uequili X, a( r) ) + ^dniu = 0 )' 7 „(x,a(r)) + c.c. 

n 

= Uequil(x,a(T)) + J ds a(s) | (^2 d n( T ~ s)in(.x,a(s)) 

+ ^2fn(T-s)q n (x,a(s)) + c.c.j 


(63) 


(64) 


It is shown in the appendix that the d ; ( t + = 0) are governed by the following equation, derived from 
equation (64): 


dj(t+ = 0 ) = [ ds a(s)( ^(^(rJJ.V^r-s)— i(i>(4)) + Y\g n (T - s) ;: ) n (x l a(s)) ) 
J — oo n n 

(65) 


where the quantities p„( t — s) obey linear integral equations involving the unknown quantities d ; ( s = 0) (/ = 
1 , . . .) and other quantities that are known such as 7 „ , u e quii, Ajrm , . . . , and D }nm . Accordingly, equation (65) 
corresponds to a Volterra integral equation of the second kind. 


To proceed further, we need to consider two cases, first a body undergoing an externally forced motion, 
and second, a body undergoing an unforced motion. In the former case, the motion history is externally 
imposed and therefore known; in the latter case, it is unknown. 
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Externally Driven Motion 


In this case, the angle of attack a( s) is known. The Volterra integral equation (65) governing dj(t+ = 
0) can be solved in principle, once the values of the equilibrium state of the velocity field u e quii . of the eigen- 
values \j(j = 1 , . . .) and of the eigenvectors 7 ; (/ = 1 , . . .) are determined from the known values of the 
angle of attack. 

Experimentalists commonly take the forcing as a periodic function of time, athough one can imagine 
different types of forcing. We shall show in the following that in the case of periodic external excitation, the 
coefficients d ; (t+ = 0) can be decomposed into a Fourier series of frequency w, the forcing frequency. Indeed, 
according to equation (65), the initial values d } (t = r)(j = 1 , . . .) have the following analytical form: 

dj(t = t) = f ds d(s)C> ; (dn(t = s)(n= \ ,a(s) ,a(r) ,t - s) j= 1,... (66) 

J — OO 


We show first that the integral constituting the right-hand side of equation (66) is periodic in time with 
same period T as the external forcing. To that end, let us make the following change of variables: 


r = t' + T 1 

s = s' + T J 


(67) 


The integral becomes, by using the periodic property of a( s ) , 

[ ds a(s)Q> } (dn(t = s)(n = \ a(s) , a(r) , t — s) 

J — OO 

I 

= J ds a(s)Q>j (dn(s + T)(n= 1, ...), or( s'), qt(t),t - s') (68) 

Therefore, according to (68) 

d n (s+T)=d n (s) (69) 

which is the property announced previously. Because the dj(t+ = 0 )(;’ = 1 , . . .) are periodic in time with 
period equal to T, they can be decomposed into a Fourier series of frequency = 2 n/T: 

dj(t + = 0) = '£d Jn e™ T (70) 

n 

Since the value of dj(t+) is governed by equation (26), it has the same property as dj(t+ = 0) of being periodic 
in time with the same frequency u>. On returning to the value of the velocity field given by equation (29), it 
is evident that u(x, a, x)) is also periodic in time with the same frequency as the excitation. Accordingly, 
the velocity field can be decomposed into a Fourier series of frequency w. This result is in agreement with 
theoretical predictions of Joseph (ref. 7) related to a body excited by a periodic forcing, in a Reynolds number 
range such that the equilibrium state at fixed a is time-invariant and stable. 

One can approximate equation (70) in the case of a slowly varying motion as 

dj(t + = 0) ~ djo + dj\e iuT + dj-i e~' UT (71) 

that is, 

dj(t+ = 0) = d } (a(r) ,a(r)) (72) 
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Elastically Mounted Body in Unforced Motion 


In the case of an elastically mounted body in unforced motion, the angle-of-attack motion history a( s ) 
is unknown. One is faced with solving equation (65) coupled with the equation of motion governing the body, 
according to a scheme of advancing step by step in angle of attack. In addition to being complicated, the above 
procedure does not give any insight into the form of dj ( t+ = 0) . It is possible, however, to get an approximate 
analytical form of dj(t+ = 0) by making use of the following reasoning. Let us first note that according to 
equations (25), (26), (54), and (65), the value of the velocity field contains terms such as 

exp{X ; (a(r))(t-r)} j- 1,... 

under an integration over time. Hence, because all the real parts of \j are negative, the effect of a distant-past 
motion history of the angle of attack on the present value of the velocity field is fading away, so that it is 
principally the recent motion history of a that determines the present value of the velocity field. In as much 
as one is concerned only with a recent-past motion history, one can consider fitting it over the recent past to 
a known motion history of the system. The latter can be provided, for instance, by experiments made under 
external periodic forcing. Under this approximation, the analytical dependence of d)-(t+ = 0) for an elastically 
mounted body will be the same as for a periodically driven body. 

For a slowly varying unforced motion, d ; (t+ = 0) obeys equation (72). As a( t) <C a(r), one can 
expand df. 

dj(U = 0) = d,{a(r),a{T)) = dj(a(r), 0) + d(r)-^- + t?(d 2 (r)) (73) 

... . doi(r) 

As a consequence of the results of this section, one can rewrite relation (62) as the following, by using relation 
(72): 

da -* 

u( x, t ) ~ Uequii( a( t) ) + J dr— a( r) , d( t) , t+ ) ( 74) 

DETERMINATION OF NORMAL FORCE 


The pressure can be derived from the velocity field by using the Navier-Stokes equations. According 
to equations (12), (20), and (21), the pressure depends on dj{U){j = 1 , . . .) as 


p( X, f +) — Pequil(. %) + Ptrans ( dj{_ t+) , ] — 1 , . . .) 


(75) 


Since the d ; (t+)(; = 1 , . . .) depend on their initial value according to equations (25) <md (26), one has 

p(M+) = PefluiiC 3-i ®( 0 ) P(rons(^/(f+ = 0)>7 = !)•••) (76) 

The value of the normal force N (normalized with respect to density p) is composed of a contribution 
from the (normalized) surface pressure and a viscous contribution arising from the skin friction at the surface. 
We derive a form for the normal force containing the two components that is applicable in the case of a 
two-dimensional flow. The restriction is imposed for simplicity only; generalization to accommodate three- 
dimensional flow is straightforward. Let S denote the body cross section, / a unit vector directed normal to 
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the axial velocity Uo , n and t unit vectors, respectively, normal and tangential to the surface element dS. Let 
C be the magnitude of the vorticity vector, determinable from the velocity field through £ =| V x u |. Then, 

N= f p(x,t)\s J.ndS + [ v <( x, t ) |s J.tdS ( 77) 

Js Js 

As shown previously, p and u (and hence Q can be separated, respectively, into an equilibrium com- 
ponent, which depends only on the instantaneous value of a, and a transient component, which depends on the 
initial values of dj{) = 1 It follows that Cn, the dimensionless form of the normal force, may be given 

the form 

C N (t) = C^ U, W» - dT^[dj[a T ( 0,U = 0] ,t + ,T ] (78) 

where the value of the function T approaches zero asf — r — ► oo, a result of the properties ofci n (t + ). Relation 
(78) was derived in previous work (ref. 8) on the basis of a functional analysis approach; in reference 8, the 
function T was called the deficiency function. The value of the normal force can be approximated as the 
following, by using equations (72): 

C N {t)~ CT\ «( 0 ) - J dr a{ r ) , d( r) , U ) ( 79) 

and by using equation (73): 

CM = - I* dT^(a(T),U) + ^ia 2 (t)) (80) 

Jo ar 

~ C% uil (a(t)) + a(t)C N ,(a(0) + t?(d 2 (t)) (81) 

According to relation (81), the value of the normal force depends only on the value of angle of attack 
a and of its derivative a at the time of observation t. In other words, only a recent-past motion history of the 
angle of attack at £ < t is taken into account through the dependence on a( 0 . The logic underlying relation 
(81) forms the basis for most mathematical models in flight dynamics. In spite of its simplicity, it has been 
capable of reproducing a great variety of the characteristics of aircraft motions (see, e.g., ref. 9). 


CONCLUSIONS 


We conclude that in the case of a time-invariant equilibrium state, derivation of the indicial response 
of the velocity field dir_„dy from the Navier-Stokes equations leads to an analytical form of the normal force 
which confirms the form of existing mathematical models. Applied to the next physical situation of a periodic 
time-varying equilibrium state, our approach will show that the same level of mathematical modeling will no 
longer hold, but rather will require considerable amendment, owing to the presence of the phase characteristic 
of the periodic equilibrium state. 
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APPENDIX 

DERIVATION OF EQUATIONS GOVERNING INITIAL CONDITIONS 


1 ,...). 


In the following, we propose to derive the equations governing the initial conditions d ; ( t+ = 0)(; = 

Let us start with the following equation derived directly from equation (64) of the main text: 
Y^dn(U = 0 ) 7 „(f,a(r)) + c.c. = J da a(s) 4 ( 7 - - s) 7 „(i, 0 ( 5 )) j 

+ ^2fn(r - s)^ n (x,a(s)) + c.c.J 


(Al) 


Multiplying equation (Al) by the adjoint vector t*( x, a( t) ) and integrating over space, one gets 

dj(t+ = 0) = [ ds a(s)(^j(x,oi(T)),^2d n (T- s)^—ln(x,a(s)) 

J-00 1 V da(s) 

+ + /n( T -s)) 7 „(i,a(a)) ) (A2) 


Equation (A2) would constitute a Volterra integral equation of the second kind if the quantities ^d ; ( r— 
a) (; = 1 , . . .) were not present. We shall show in the following that, indeed, the quantities g ; (j = 1 , . . .) 
defined as 

9 j = ( T - s) + fj ( t - s) ( A3) 

do not contain any element such as ^d ; (r — a)(; = 1 , . . .). 

Let us first write the following equations governing fj ( t — s), derived from equations (54): 

/;( T - a) = - JLdj.e x * T - a) + (vj + j-dj - Cj H( r - a) ) e x ' (r -’> 

+ J2 A i™ [ T ~ 3 dsie^ T ~*~ tl ^ (dn(si) f m (s\) + d m (si)/ n (si)) + ... 

+ £ Djnm f T ~ S cfeie X/(T- *“ 5,) (3»<*l>7mC*i) + 3»(*l)>„(«i)) (A4) 


n^m 


where the quantities + j^d/j do not contain any element ^d ; (r- s)(/ = 1,...). Indeed, one can show 
by using equation (51) that 

3 d 

V] + -^dj = -( Tj(x,a(s)), -^—^ii equ ii(x,a(s)) 


da( a) 


+ ?(« ( *= 0) ^( ? '“ W> + 3 ”8^)’") > <A5) 
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The quantities = 1 , . . .) can be obtained from equation (24) as 


JL d = 

da } 


(T ~ s) 'L d ’ (s)+ ^ Ajim Io s<i31 eA/(T ' s_,,) (^ dn(5i) dTn(5 1) 

■d n( S,).^d m ( S l ) ) + ...+ X)0 ; n ro ^ T_S ( i Sl *> (J^d n(si).d m (si) 

+ Qj 




(A6) 


with Q ; (y = 1 , . . .) defined as 

d 


Qj = dj(s).j-e x * T - a) + <ki (A Jim e X i (T - a - ai) ) d»(«i)dm(ai) 


n,r7i 

••. + £ f S ^(Dj^e^- 8 - 8 ^) InisOlrnCsQ j= 1,... 


W J ° 


dcx 


(A7) 


that is, the quantities Q ; are independent of -g^djij = 1, ...). 
Adding equation (A4) to (A6), one gets 


9, = -'>«<♦>) 

+ £ A/nm f dsi e x ' (T ~ a ~ Sl) (dn(si)g m (si) + d m (si)g n (si)) + ... 

n,m 

+ E D ;»m / T S dsi e^ (T -*- Sl) (^(sOp^Sl) + dm(8i)p B (ai)) . ; = 1,... (A8) 

n ,m 70 

The equations governing gj(j = 1 , . . .) do not contain any element -s^dji) = 1 as announced previously. 
Thus, equation (A2) does in fact constitute a Volterra integral equation of the second kind, as stated in the text. 


18 



REFERENCES 


1. Tobak, M.; Chapman, G. T.; and Unal, A.: Modeling Aerodynamic Discontinuities and Onset of Chaos 

in Flight Dynamical Systems. Ann. Tel6commun„ tome 42, no. 5-6, 1987, pp. 300-314. Also NASA 
TM-89420, 1986. 

2. Sarpkaya, T.: Vortex-Induced Oscillations. J. Appl. Mech., vol. 46, no. 2, 1979, pp. 241-258. 

3. Bearman, P. W.: Vortex Shedding from Oscillating Bluff Bodies. Ann. Rev. Fluid Mech., vol. 16, 1984, 

pp. 195-222. 

4. Yates, L.; Unal, A.; Szady, M.; and Chapman, G. T.: One-Degree-of-Freedom Motion Induced by Mod- 

eled Vortex Shedding. NASA TM-101038, 1989. 

5. Joseph, D. D.: Stability of Fluid Motions I. Springer Tracts in Natural Philosophy, vol. 27, 1976, pp. 25-62. 

6. Fliess, M.; Lamnahbi, M.; andLamnahbi-Lagarrigue,F.: An Algebraic Approach to Nonlinear Functional 

Expansions. IEEE CAS, vol. 30, no. 8, 1983, pp. 554-570. 

7. Joseph, D. D.: Remarks about Bifurcation and Stability of Quasi-Periodic Solutions which Bifurcate from 

Periodic Solution of the Navier-Stokes Equations. Lecture Notes in Mathematics 322. 1. Stakgold, 
D. D. Joseph, and D. H. Sattinger, eds., Springer 1973, pp. 130-158. 

8. Tobak, M.; and Pearson, W. E.: A Study of Nonlinear Longitudinal Dynamic Stability. NASA TR R-209, 

1964. 

9. Hui, W. H.; and Tobak, M.: Bifurcation Theory Applied to Aircraft Motions. AGARD Symposium on 

Unsteady Aerodynamics: Fundamentals and Applications to Aircraft Dynamics, AGARD CP-386, 
1985, pp. 26-1 to 26-13. 


19 



Report Documentation Page 

Space Admrtsfretan 

1. Report No. 

NASA TM- 102856 

2. Government Accession No. 

3. Recipient's Catalog No. 

4. Title and Subtitle 

Indicial Response Approach Derived from Navier-Stokes 
Equations 

Part 1 - Time-Invariant Equilibrium State 

5. Report Date 

October 1990 

6. Performing Organization Code 

7. Author(s) 

K. V. Truong and M. Tobak 

8. Performing Organization Report No. 

A-90258 

10. Work Unit No. 

505-61-71 

9. Performing Organization Name and Address 

Ames Research Center 
Moffett Field, CA 94035-1000 

11. Contract or Grant No. 

1 3. Type of Report and Period Covered 

Technical Memorandum 

1 2. Sponsoring Ayency Name and Address 

National Aeronautics and Space Administration 
Washington, DC 20546-0001 

14. Sponsoring Agency Code 


15. Supplementary Notes 

Point of Contact: 


M. Tobak, Ames Research Center, MS 227-4, Moffett Field, CA 94035-1000 
(415) 604-5855 orFTS 464-5855 


/ 


/* 




r 4 


OtS 1 




1 6. Abstract 


c>fc« 


the indicial response approach in a form appropriate to the study 
of vortex-induced oscillations phenomena. We-demenstrate-that an appropriate form tor the indicial 
response of the velocity field may be derived directly from the Navier-Stokes equations. Fet-conveaiencer- 
4he-stedy^s-divided-interthree^arts i -In-Rart4rwe-demonstraKr, on the basis of the Navier-Stokes equations, 
‘j**how a form of the velocity response to an arbitrary motion may be determined. To establish its connection 
with out previous work, the new approach is applied first to the simple situation wherein the indicial 
response has a time-invariant equilibrium state. Results for the aerodynamic response to an arbitrary 
motion are shown to confirm the form th at - we - h ave obtained previously. 


1 7. Key Words (Suggested by Authors)) 

Unsteady aerodynamics 
Nonlinear flight dynamics 
Aerodynamics 

Nonlinear dynamical systems 


18. Distribution Statement 

Unclassified-Unlimited 
Subject Category - 02 


19. Security Classif. (of this report) 

20. Security Classif. (of this page) 

21. No. of Pages 

22. Price 

Unclassified 

Unclassified 

22 

A02 


NASA FORM 1626 OCT 86 


For sale by the National Technical Information Service, Springfield, Virginia 22161 




